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J> | Recent developments of the small x CCFM evolution are described, in- 

^^O . eluding improvements of the splitting function. The resulting unintegrated 

gluon densities are used for predictions of hadronic final state measure- 
ments like jet production at HERA and heavy quark production at the 
Tevatron. 
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C*~) ■ This paper is dedicated to the memory of Jan Kwiccinski, who passed 

away so early. 

Oh' 

1. The problem with collinear factorization 

Ph! __ 

D , Inclusive cross section measurements [2121; ^ut a ^ so man y hadronic fi- 

^ | nal state cross sections 011103 are in general well described in QCD. The 

cross sections are factorized into integrated parton density functions (PDFs) 
which are evolved from a starting scale fio to the factorization scale \x / using 
the DGLAP jfHEllElini evolution equations and convoluted with the process 
dependent coefficient functions (matrix elements). This approach is often 
called collinear factorization. 

However, regions in phase space exist where the above mentioned ap- 
proach using coefficient functions evaluated at fixed order in the strong 
coupling constant, a s , is not sufficient to describe the measured cross sec- 
tions. Such measurements are e.g. forward jet production and the azimuthal 
de-correlation of jets at HERA, but also heavy quark production at the 
Tevatron. 

Processes sensitive to hard multi-parton emissions are generally not cov- 
ered fully in DGLAP. The BFKL [TOUnUH] approach, also called the k t - 
factorization [131114] or the semi-hard approach [151116] . is assumed to be 
better suited for such a scenario, since partons in the initial state cascade 
can have any kinematically allowed transverse momentum, in contrast to 
DGLAP. The CCFM HTHTHUmHOI approach and its reformulation in the 
Linked Dipole Chain Model LDC [2TU22] attempt to cover both the DGLAP 
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and BFKL regions by considering color coherence effects. A general intro- 
duction to small-x physics and the evolution equations can be found in [23] ■ 

2. CCFM - LDC evolution and unintegrated parton densities 

In the CCFM evolution equation angular ordering of emissions is intro- 
duced to correctly treat gluon coherence effects. In the limit of asymptotic 
energies, it is almost equivalent to BFKL [2Hl25U26| . but also similar to the 
DGLAP evolution for large x and high Q 2 . The cross section is fe-factorized 
into an off-shell matrix element convoluted with an unintegrated parton 
density (uPDF), which now also contains a dependence on the maximum 
angle H allowed in emissions. This maximum allowed angle H is defined by 
the hard scattering quark box, producing the (heavy) quark pair and also 
defines the scale for which parton emissions are factorized into the uPDF. 

The original CCFM splitting function is given by: 

a s {p t ) a s {k t ) f 1 dz' [ k u dq 2 _ 

Pg = -, 1 A ns with logA ns = -/ — / ~^-a s (1) 

1-Z Z J Zl Z' J( z 'q t )2 q 2 

with a s = The phase space region of angular ordered emissions is 
obtained from: 

Zi-iqi-i < q~i with % = — (2) 

1 Zj 

Here Z{ = Xi/xi-i is the ratio of the energy fractions in the branching 
(i — 1) — > i and pu is the transverse momentum of the emitted gluon i. 
The transverse momentum of the propagating gluon is given by ku- It is 
interesting to note, that the angular ordering constraint, as given by eq. ([%)). 
reduces to ordering in transverse momenta pt for large z, whereas for z — > 0, 
the transverse momenta are free to perform a so-called random walk. 

The LDC model is a reformulation of CCFM, where the separation be- 
tween the initial- and final-state emissions is redefined and, in addition to 
the angular ordering in eq . (|2"]l . the gluons emitted in the initial state are 
required to have 

p ±i > min(Xu, fey-i). (3) 

This constraint means that the pt of the emitted gluon is always close to 
the highest scale in the splitting and the argument in a s is naturally taken 
to be pt- 

While formally equivalent to the DLLA accuracy, it is important to note 
that the sets of chains of initial-state splittings summed over, are different 
in LDC and CCFM. One chain in LDC corresponds to a whole set of chains 
in CCFM, and predictions for non-inclusive observables of the evolution 
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will only be comparable if the correct final-state emissions are added in the 
regions relevant for the formalism used. 

The Cascade [271125] program is a direct implementation of the CCFM 
approach into a hadron level Monte Carlo generator using a backward evolu- 
tion approach for the initial state parton shower. The Ldc |2~9TI30) program 
is a realization of the LDC model in terms of a hadron level Monte Carlo 
generator. 



2.1. Unintegrated gluon densities 

The original CCFM splitting function given in eq.Q includes only the 
singular terms as well as a simplified treatment of the scale in a s . Due 
to the angular ordering a kind of random walk in the propagator gluon kt 
can be performed, and therefore care has to be taken for small values of 
kf. For values of kt < k1 ut the non-perturbative region is entered, which is 
avoided in a strictly pt-ordered evolution (DGLAP). In CCFM, for kt < kt ut , 
a s {kt) and also the splitting probability could become un-physically large. 
However, this is the region of high parton densities, where saturation effects 
could appear. A practical treatment here is to freeze a s and the gluon 
density for values k t < kt ut . 

Three new sets of uPDFs were determined [H^: J2003 set 1 with the 
splitting function given in eq. Q , J2003 set 2 including also the non-singular 
terms in P g and J2003 set 3 using pt as the argument in a s . For all sets 
the input parameters were fitted to describe the structure function F2 as 
measured at HI [321133] and ZEUS |341l35j in the range of x < 10 -2 and 
Q 2 > 5 GeV 2 . Using 248 data points a x 2 /ndf = 1.29, 1.18, 1.83 for J2003 
set 1,2,3 , respectively, is obtained. A comparison of the different sets of 
CCFM uPDFs is shown in Fig. [Q 

Also the LDC model describes F% satisfactorily well, but the correspond- 
ing unintegrated gluon densities are somewhat different, as also quarks can 
be included in the evolution. In Fig. |2] three different unintegrated gluon 
densities for the LDC approach are presented. The standard set refers to 
the full LDC including quarks in the evolution and the full gluon splitting 
function, whereas for the gluonic set and the leading set only gluon evolu- 
tion is considered with only singular terms in the splitting function for the 
latter. All three alternatives have been individually fitted to F% in the region 
x < 0.3, Q 2 > 1.5 GeV 2 for standard and x < 0.013 and Q 2 > 3.5 GeV 2 for 
gluonic and leading. In LDC there is only one relevant infrared cutoff, kto, 
which limits the pt of emitted gluons. 
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Fig. 1. Comparison of the different sets of unintegrated gluon densities obtained 
from the CCFM evolution as a function of x for different values of kt at a scale of 
q = 10 GeV. 



3. Comparison with hadronic final state data 

A comparison of measurements of hadronic final state properties, like jet 
or heavy quark cross sections, with theoretical predictions requires Monte 
Carlo event generators, which also allow to apply the hadronization step. In 
the following sections a few examples are presented where data are compared 
to predictions obtained with Cascade based on the CCFM unintegrated 
gluon densities, described in section 12.11 



3.1. Jet cross section at HERA 

The azimuthal correlation of dijets at HERA is sensitive to the trans- 
verse momentum of the partons incoming to the hard scattering process and 
therefore sensitive to the details of the unintegrated gluon density. This was 
studied in a measurement [HOI °f the cross section for dijet production with 
Et > 5(7) GeV in the range 1 < r/i a b < 0.5 in deep-inelastic scattering 
(1CT 4 < x < 1(T 2 , 5 < Q 2 < 100 GeV 2 ). In LO collinear factorization, 
dijets at small XBj are produced essentially by — > qq, with the gluon 
collinear to the incoming proton. Therefore the qq pair is produced back- 
to-back in the plane transverse to the 7*p direction. From NLO (C(a 2 )) 
on, significant deviations from the back-to-back scenario can be expected. 
In the /^-factorization approach the transverse momentum of the incom- 
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Fig. 2. Comparison of the different sets of unintegrated gluon densities obtained 
within LDC at scale of q — 10 GeV. standard refers to the full LDC including quarks 
in the evolution and the full gluon splitting function. For gluonic and leading only 
gluon evolution is considered with only singular terms in the splitting function for 
the latter. Also shown is the J2003 set 1 for comparison (divided by it). 



ing gluon, described by the unintegrated gluon density, is taken explicitely 
into account, resulting in deviations from a pure back-to-back configuration. 
The azimuthal de-correlation, as suggested in Ref. [37], can be measured: 

s= ffiWA^ow 0<a<180 . 

Jo 180 N 2 - jet (Ap,x,Q 2 )dAp 

In the measurement shown in Fig- El a = 120° has been chosen. The data are 
compared to predictions from Cascade using J2003 set 1-3. Also shown 
for comparison is the NLO-dijet J3S calculation of the collinear approach. 
One clearly sees, that a fixed order NLO-dijet calculation is not sufficient, 
whereas J2003 set 2 gives a good description of the data. However, the 
variable S is sensitive to the details of the unintegrated gluon distribution, 
as can be seen from the comparison with J2003 set 1 and set 3. 

A measurement, aiming to observe deviations from the collinear DGLAP 
approach, is the production of jets in the forward (proton) region. The phase 
space is restricted to a region of Q 2 > 5 GeV 2 and E^ 1 > 3.5 GeV in 
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Fig. 3. The S distribution for dijet 
events (SHj- The prediction from Cas- 
cade are shown together with the 
NLO-dijet prediction in the collinear 
factorization approach. 
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Fig. 4. The cross section for forward jet 
production as measured by HI |39j as a 
function of x Bj ■ The prediction are the 
same as in Fig. El 



the forward region of 1.7 < rjj e t < 2.8 with the additional requirement of 
0.5 < E\ j et /Q 2 < 2, a region where the contribution from the evolution in 
Q 2 is small. The cross section for forward jet production has been measured 
by HI as a function of XBj, shown in Fig. together with predictions 
from Cascade. Also the NLO-dijet |M8| prediction in the collinear approach 
is shown. The fixed NLO-dijet calculation falls below the measurement, 
whereas the ^-factorization approach supplemented with CCFM evolution 
gives a reasonable description of the data. 



3.2. Heavy Quark Production at the Tevatron 

The differential cross section as a function of the transverse momentum 
of D-mesons has been measured in pp collisions at = 1.96 TeV by the 
CDF collaboration |40| . They find the measured cross section to be larger 
than the NLO predictions in the collinear factorization approach by about 
100 % at low pt and 50 % at high p%. In Fig. [5] the measurement is shown 
together with the predictions obtained in /^-factorization using Cascade 
with the CCFM unintegrated gluon density described above. For J2003 set 
1 and set 3 good agreement for all measured charmed mesons is observed. 

The cross section for bb production in pp collision at ^/s = 1800 GeV 
has been compared with the prediction of Cascade based on the CCFM 
gluon densities in jlj. It is interesting to note that J 2003 set 2 which 
describes best the data of Fig. |31 predicts a lower cross section for heavy 
quark production. 
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Fig. 5. Differential cross section for D meson production as measured by CDF |^0] 
as a function of the transverse momentum compared to predictions from Cascade. 



4. Conclusion 

It has been shown, that fcj-factorization and the CCFM evolution of the 
gluon density is a powerful tool for the description of hadronic final state 
measurements. Improvements of the CCFM splitting function have been 
discussed. 

Jet measurements at HERA, but also measurements of charm and bot- 
tom production at the Tevatron can be reasonably well described, whereas 
calculations performed in the collinear approach even in NLO have dif- 
ficulties to describe the data. This shows the advantage of applying kt- 
factorization to estimate higher order contributions to the cross section but 
also the importance of a detailed understanding of the parton evolution 
process. 
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